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ABSTRACT 

'a\ '. 
ov 

Intrinsically cold particle dark matter inevitably creates halos with sharp 
discontinuities in projected surface density caused by the projection of fold 
catastrophes onto the sky. In principle, these imperfections can be detected and 
■ measured with gravitational lensing through discontinuities in image magnification 

, and image structure. Lens solutions are discussed for the most common universal 

CN , classes of discontinuities. Edges caused by cold particles such as condensed axions 

\ and thermal WIMPs are very sharp, respectively about 10 -12 and 1CP 7 of the halo 

scale. Their structure can be resolved by stellar and quasi-stellar sources which show 
' sudden changes in brightness or even sudden disappearances (sometimes within hours) 

CN ■ as edges are crossed. Images of extended objects such as edge-on galaxies or jets can 

show sudden bends at an edge, or stretched, parity-inverted reflection symmetry about 
a sharp line. Observational strategies and prospects are briefly discussed. 

G\ ■ 

^ , 1. Singularities in Projected Density 

£ : 

Particle dark matter naturally forms thin structures in phase space QSikivie 1999|) . In some 
cases, such as axions forming by Bose condensation, the peculiar velocities are exceptionally small, 
• i-H ! but even for thermally produced particles, adiabatic cooling creates a cool Hubble flow which 

resembles a 3D sheet (v = Hr) in six-dimensional phase space. The fine-grained phase structure 
remains thin even after nonlinear collapse destroys the ordered Hubble flow. For systems which 
are not too dynamically old or too violently relaxed, the fine structure retains information about 
the nature of the dark matter and its dynamical history. 

The fine-grained structure remains most conspicuous at the locations of caustics in the flow of 



matter. The most famous examples in cosmology are the Zel'dovich "pancakes" ( Zeldovich 1970 ). 
Similar caustics occur in general nonlinear collapses in smooth potentials: for example, shells in 
spherical similarity solutions (Fillmore & Goldreich 1984, Bertschinger 1985| ), or ringlike tubes in 



collapses having less symmetry ( Sikivie 1999| ). Singularities are not an artifact of the symmetry: 



fold catastrophes (2D sheets) joining in cusp catastrophes (ID lines), which themselves join in 
various pointlike catastrophes, are a generic feature of any smooth 3D phase sheet mapping onto 
3D space ( [Arnold 1984 , Tremaine 1998| ). Physically these structures are caused by turning points 



or focus points of orbits. They create formal singularities in density (sheets, lines and points) 
with universal density profiles. These features may be observable through their dynamical effects 
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( Sikivic 199S| ), though with insufficient resolution to reveal the fine structure of the dark matter 



phase sheet. 

Here I discuss the projection of the phase sheet onto the sky and how discontinuities in surface 
density can manifest themselves in gravitational lensing. Lensing offers the advantage of very 
small background sources which can serve as microscopic probes of small scale phase structures. 
Observations of such events can yield physical information about the nature of the dark matter 
and its formation mechanism- for example, whether axions formed by condensation or string 
radiation. 

The most common generic true spatial catastrophes of cold dark matter — a 3D phase sheet 
moving in a 3D subspace of 6D phase space — are "folds", 2D structures of codimension and 
corank 1 ( [Tremaine 1998 ). They terminate on linelike cusps which in turn terminate on pointlike 



swallowtail, elliptic umbilic and hyperbolic umbilic catastrophes. Folds form sheets with a formally 
divergent density, with a universal density profile: on one side p oc x -1 / 2 , abruptly falling to p = 
on the other side.Q However, when projected onto the sky the catastrophic divergence does not 
survive. Integrating the fold profile in projection produces instead a a simple jump in column 
density, with a sharp edge, along the line where a fold sheet is convex and tangent to the line of 
sight. 

Consider a (generic) 2D fold of a 3D phase sheet in 3D space. Let the z axis lie along the line 
of sight and tangent to the fold, the y axis in the plane of the sky and tangent to the projection 
of the fold at x = 0. The sheet lies at x s (z) = z 2 , so the 3D density to the right of the sheet is 
p = (x — z 2 )~ x l 2 . The projected density is 

f=+^ j ir forx>0 

P(X) = L-V, P{Z )dZ = \ for x < 0, (1) 

a discontinuous jump where the projected fold is tangent to the line of sight. The projected 
edges form lines on the sky with coherence scale £, along which the amplitude of the discontinuity 
slowly varies. The amplitude of the jump (which depends on the curvature of the fold in the 
invisible dimensions), and its coherence scale on the sky, depend on the dynamical substructure 
and history of the halo. Roughly speaking, a "smooth" halo which is N orbital times old has 
about N wraps and folds of its phase sheet, and about iV projected edges. A typical discontinuity 
has an amplitude <5E ~ A^ _1 S where the mean halo surface density is S. 

In reality the catastrophes are concentrated in the central, dynamically old parts of a halo 
and display a chaotic structure reflecting the dynanical history of the system; only the younger 
folds, representing particles on their first or second orbits, show coherence on the scale of the halo 



Note that these structures actually exist in the fine-grained halo density distribution. In principle, a detailed 
time history of the WIMP flux on Earth would show spikes of divergent flux every few million years. The integrated 
flux in these spikes however is small so in a real fossil record (such as tracks in mica) they would not stand out clearly 
against the background. 
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itself. Although the statistical properties of the edges in realistic models have not been studied 
in detail, they are apparent in 2D projections of dynamically young systems in high-resolution 
N-body simulations (J. Stadel, private communication). Edges disappear if the matter is too hot 
or if the experiment is too coarse-grained compared to the scale of the dynamical wraps. 

At the finest resolution, the sharpness of the edge discontinuities is in principle limited only 
by the intrinsic width of the dark matter phase sheet. This scale is typically much smaller than 
the angular scale t of the dynamical folding, in which case the lensing effects can be treated as a 
one-dimensional problem. 



2. Lensing by Projected Folds 
2.1. Lenses with ID Translational Symmetry 



We follow the notation of Blandford and Narayan (1986), describing relativistic bending 
in the thin-screen approximation by a two-dimensional potential ip on the plane of the sky 
(x,y). We study line discontinuities on the y axis and because of the translational symmetry the 
two-dimensional Poisson equation reduces to one dimension, 



dx\ 



(2) 



The critical surface density S c = c 2 (47rG-D) _1 = 0.35g cm~ 2 (D/lGpc) _1 is determined by the 
combination D = DdDd s /D s of the angular diameter distances to the lens Dd and the source D s , 
and between them Dd s ] at this density a uniform sheet of matter would just focus radiation from 
the source at the observer. The positions of the images xj are determined by the lens equation 

dtp(xi) 



xi = x s + 



dxi 



(3) 



where x s is the direction of the source. The image distortions and magnifications are determined 
by the extrinsic curvature tensor 



K, 



de 



Si 



dd 



K + H 









(4) 



where the principal axes for lensing are the same as those set by the discontinuity, 9\ = x and 
02 = y- Here k = 1 — (S/S c ) denotes the expansion, and \x denotes the shear associated with image 
distortion. By symmetry in this situation, the yy term of K is unity, determining a particular 
relation between expansion and shear: n — fj, = 1, and hence \i = — (S/S c ). The magnification of 
an image is given by M = \k 2 — /U 2 )" 1 and hence 



M 



S 



(5) 



images are stretched in the x direction by this factor. 
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2.2. Folded 2D phase sheet 



Some dissipative structures (such as cold disks) resemble a 2D phase sheet, whose projection 
onto the sky can create high density features resembling a fold catastrophe in projection (for 
example, stellar shells around elliptical galaxies; |Quinn 1984| , |Hernquist & Spergel 1992]) . In 
reality these sheets are not as thin as intrinsically 3D dark matter phase structures although they 
may produce recognizable pseudo-catastrophe events. 



In the plane of the sky consider a fold on the y axis at x = 0. The universal surface density 
profile is 

E(x) _ f x" 1 / 2 for x > 
E " for x < 



(6) 



We have used the scaling freedom to set £ = S c at x = 1. The solution of the lens equation yields 
the mapping between source and image positions, 

(2± V4 + x 5 ) 2 for xi > 
xs for xi < 



X/ 



(7) 



Three images in a line are produced for x,5 in the vicinity of the fold, —4 < x$ < 0. The 
magnifications of the images at xj > are 



M 4 



II - 2x 



-1/2,-1 



1 ± 



1 + T 



-1/2 



(8) 



These formulae describe a universal family of solutions characterized by one parameter. However, 
the approximation breaks down on a scale determined by the temperature of the 2D phase sheet; 
while it might serve as a good description of some baryonic lenses, we cannot use this as a probe 
of dark matter. 



2.3. Sharp Edge from Folded 3D Phase Sheet 

The more interesting case is the projection of generic catastrophes caused by folds of the 
very cold 3D CDM phase sheet, which create discontinuities with extremely sharp edges, and a 
projected profile 

£(x) J (Tq for x > 



S c for x < 0. 



(9) 



The lensing solution is simply 




, .. „ v .- - 2<7 ) for xi > 
xi = { (10) 

for xi < U 



There are no multiple images for 2<7o < 1. The mapping from S to I receives a stretch in the x 
direction to the right of the edge. There is an abrupt change in magnification from one side of the 
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edge to the other, from M = 1 to M = |1 — 2<tq| , which would appear as a sudden change in the 
brightness of a small moving source. The discontinuity can also be observed in image structure if 
there is an extended background source (such as a compact jet or edge-on disk) which lies in a 
straight line: at the location of the edge the image displays a sudden bend. 

If 2(Jo > 1, two images form if x$ < (at xi > and xj < 0) and no images form if x$ > 0. 
A source moving towards the edge from x < appears as two images approaching each other with 
a uniform magnification (one image with M = 1, the other with M = |1 - 2cr |~ 1 , both moving 
with angular velocity oc M) until it reaches the edge, when the images meet and both suddenly 
disappear. The image of an extended object at xs < appears reflected across the line, stretched 
and parity-reversed; that is, the surface brightness at (x,y) is duplicated at (—Mx,y). 

These very simple solutions only apply over angles much smaller than £; particularly in the 
supercritical case there are nontrivial macroimages with a larger-scale structure not described 
here. However, these solutions accurately illustrate the unique effects near the edges. 



3. Observing the Innermost Phase Structures 

Physically, the edge discontinuity is the most interesting case since the width of the edge 
reveals the initial dark matter velocity dispersion. An extreme example is condensed axions, with 
dispersion at the redshift of collapse z co u ( [Bikivie 1999 ) 

5v a /c «3x 10- 17 (m a /10- 5 eV)-\l + z coU ). (11) 

A more typical case is thermally produced WIMPs which have 

5v w /c « lO- n {m w /GeV)V 2 (l + z coll ). (12) 

This dispersion creates a range of energies at a given initial position, translating after one orbit 
to a dispersion of turnaround radii for the first fold and smearing the projected edge to a fraction 
of the halo size ~ 8v/v v i r where v v i r (~ 10 _3 c) is the halo dispersion. This ratio is very small — 
adopting z co u ~ 3 to 10 typical of CDM models, it is about 10~ 12 for axions and 10 -7 for WIMPS. 
If the edge results from a fold of angular size £, its angular width is about 6x ~ £5v/v v i r . The 
edges are therefore very sharp, ranging from a few milliarcsec for WIMP edges in our local halo 
(£ « 0.1, 6x ps 10" 8 )) to a few picoarcsec for axion edges in clusters or halos at high redshift 
(£ m 10" 4 , 5x w 10" 16 . 

Resolving the structure in these extraordinarily sharp edges — the interesting goal if we want 
to use them to probe dark matter — requires small (and hence faint) background sources. For 
halo-scale folds (~ 10 23 cm), normal solar-type stars (~ 10 11 cm) are about the right size to 
resolve the axion edge. Bright sources such as quasars and supernovae (~ 10 15 cm) are too large 
to resolve the axion edge but have the advantage that they are visible at high redshift, and can 
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still resolve edges produced by thermal WIMPs. In principle at least, realistic background sources 
can probe the initial dark matter distribution function. 

How practical is it to use these occurrences as a probe? One type of experiment would be 



MACHO-style (Alcock et aL 1997, Paczynski 1996): monitor a large number of background 



sources (such as a galaxy of stars) and watch for sources crossing edges in a foreground halo. 
Clearly, the light curves caused by edge transits are very different from the "standard" MACHO 
events of point-mass microlensing ( Paczynski 1986j ). The number of edges across the face of a halo 



is about the number of orbits since it formed, say N ~ 10 to 100. Each background star traverses 
the angular scale of the halo in a dynamical time thalo ~ 10 s yr so monitoring 10 8 stars for a year 
will show 10 to 100 of them crossing edges. The timescale for the brightness to change (assuming 
the source resolves the edge) is the transit time of the edge « th a i (5v /v V i r ), about an hour for the 
axion edge and years for thermal WIMPs. 

Unfortunately the amplitude oo ~ S/A r of the jumps is less than 10~ 6 for Local Group 
experiments, which is not observable given limitations of photometric accuracy and the stability 
of the background sources. To find large-amplitude or "appearing/disappearing" events, sources 
must be monitored at high redshift where do in typical halos can approach unity. 

This requires surveys which are simultaneously larger, deeper, and better sampled in time 
than is possible with current telescopes. If the goal is to resolve the axion edge, one must detect a 
solar-size star in an hour at the Hubble length, which takes a 1000m telescope (fortunately, this is 
also large enough to resolve a galaxy this far away into stars). With a 100m telescope — which is 
at least being contemplated — the experiment could be done with giant stars as targets, the events 
lasting of the order of days. 

In the short term it is more realistic to study these effects via imaging. Discontinuities can 
be found by sharp bends in images of extended linear objects, such as jets or edge-on galaxies. 
Imaging of jets with mmVLBI (with resolution potentially as high as Sx ~ 10~ 10 ) or VLBA 
(8x ~ 10 -8 ) can verify the presence of sharp edges and test for the cold character of the primordial 
phase sheet. Although the accessible angular scales are not nearly small enough to probe axions, 
they can reach the WIMP scale and test for warmer forms of dark matter. 

Most interesting are high amplitude supercritical edges (2o"o > 1) in dense halos and clusters 
at high redshift with S > S c . They create double images of extended objects, reflected about a 
sharp line and stretched perpendicular to the line. A jet would not only appear to bend sharply, 
it would appear the same (but mirror-reversed) on the two sides of the bend except for a constant 
stretch factor in x (and no stretch in y). The sharp character of the edge is apparent if the 
resolution is better than the coherence scale £ of the folds, which may be almost as large as the 
macrolens itself in dynamically young systems. A fresh fold in a young cluster, many arcseconds 
across, could reflect a whole galaxy image. Edge lensing might even appear in high-resolution 
images of galaxies and quasars in known macrolenses (e.g. Kochanek et al. 1998| ). Note that 



lensing conserves surface brightness so there must be image structure on scales below £ for the 
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edge to be visible. 

Even though it may be possible to resolve the predicted scale of the thermal WIMP candidates, 
we will have to get lucky to find a source this small lying right behind an edge. One can imagine 
searching for jets behind the small fraction of the sky (about 10 -3 ) with supercritical surface 
density in a foreground halo. Suppose we are looking at jets with fraction ~ Xj/l will 

happen to transect a young edge in the halo. This means that about a fraction 10~ 3 Xj/£ of all jets 
will be suitably placed; if we imagine that there are 10 8 jets on the sky (judging from the total 
number of galaxies and their duty cycle as AGN), we ought to be able to find 10 5 Xj/£ of them 
(times whatever aspect ratio they have) which probe edges on the scale 5x ~ Xj — or about 10 -2 
for WIMPs. So, it is realistic to expect to probe well below the halo dispersion for the coldness of 
the dark matter, but probably not reasonable to test models of WIMPS. 

The ultimate goal of using this effect to probe the primordial phase structure of particle CDM 
is therefore still some way off. However, relatively coarse resolution still offers a test of scenarios 



designed to smooth halo substructure or soften halo cores (e.g., Ghigna et al. 1998, Moore et 



al. 1997): for example, reducing small-scale fluctuation power would increase the fold coherence 



scale, and raising Sv would increase the thermal smearing scale, ultimately eliminating the edges 
altogether. Observation of a sharp dark matter edge on sub-halo scales could rule out hot dark 
matter candidates. Although the theoretical predictions need to be sharpened, it seems likely that 
some of these tests can be performed using current technology, even using background galaxies as 
sources, with a sufficiently large survey. 
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Fig. 1. — Illustration of lensing by a fold catastrophe. At top, a halo has an edge (one of many) 
with angular size i. A surface density cut across the edge is shown at two different locations, with 
different amplitudes. The graphical solution of the lens equation is shown for supercritical (left) 
and subcritical (right) cases, showing the mapping of the source position to image position. The 
appearance on the sky is shown at the bottom;a moving source is shown at four times, showing 
its disappearance in the supercritical case; an arrow shows a reflected/stretched and bent image 
respectively in the two cases. 



